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An isotropic passive scalar field T advected by a rapidly-varying velocity field is studied. The tail 
of the probability distribution P(0, r) for the difference 6 in T across an inertial-range distance 
r is found to be Gaussian. Scaling exponents of moments of 9 increase as y/n or faster at large 
order n, if a mean dissipation conditioned on 9 is a nondecreasing function of \9\. The P(9,r) 
computed numerically under the so-called linear ansatz is found to be realizable. Some classes of 
gentle modifications of the linear ansatz are not realizable. 



In the past several years, a number of papers have dealt 
with the theory of a statistically isotropic passive scalar 
field advected by an incompressible velocity field that 
varies very rapidly in time (white velocity field) [|ij-p6||. 
Some exact statistical relations can be written for this 
problem, which leads to the hope that, for the first time, 
anomalous scaling exponents for the inertial range can 
be derived analytically in a turbulence problem. 

Nevertheless, progress has been limited. Closed equa- 
tions exist for each order of n-point moments of the scalar 
field H , but above 2nd order the number of independent 
variables is daunting. Perturbation analysis has yielded 
predictions of finite-order scaling exponents at two re- 
mote borders of the problem: infinite space dimension- 
ality and infinitely nonlocal interaction of spatial 
scales of the velocity and scalar fields, with diffusivity 
exponent £ -> in (10) below §@. Neither regime is 
accessible to test by direct numerical simulation (dns). 

A physically-motivated "linear ansatz" for dissipation 
terms in the equation of motion jl 10 predicts all expo- 
nents for the full domain of d and £. It has quantitative 
support from dns |l0|,^7| and from some experiments Q . 
The ansatz has not been derived analytically. The pre- 
dictions in the limit £ — + are implausible and in conflict 
with the perturbation results. It is unclear whether the 
ansatz has any domain of exact validity. 

In the present paper it is deduced that, contrary to 
common expectation, the tail of the probability distri- 
bution function (pdf ) of spatial scalar-field differences is 
Gaussian in the inertial range, rather than exponential or 
stretched-exponential. This leads to a relation between 
the asymptotic behavior of scaling exponents and that 
of a conditional mean of dissipation. The linear ansatz 
is found to yield a realizable (everywhere positive) pdf, 
while some gentle deviations from it do not. A general 
apparatus is formulated for relating scaling exponents to 
the pdf shape in the inertial range. It is used to obtain 
the linear-ansatz pdf explicitly. 

The passive scalar field T(x, t) obeys 



9 , s „ 
^ + u(x,t).V 



where k is molecular diffusivity. A statistically steady 
state may be maintained by adding an appropriate source 
term to (Q) 

The structure functions are defined by S q (r) = 
(\9(r)\i), where 0(r) denotes T(x + r,i) - T(x,t) and 
() denotes ensemble average over homogeneous, isotropic 
statistics. Suppose that there is power-law inertial-range 
scaling of the structure functions S q (r), 



S q (r) ex (r/L) 



(2) 



where I > r > ^, L is a large length scale, the or- 
der of initial or injection scales, and id is a dissipation 
scale. A realizable (everywhere positive) pdf P(9, r) for 
9(r) requires that d^q/dq be a non-increasing function of 
q (Holder inequalities) . This does not exhaust the realiz- 
ability conditions on (^ q . Necessary and sufficient condi- 
tions for an everywhere non-negative P{9, r), nonzero at 
an infinite set of 9 values, are [p8[ 

det[S , i+J (r)] iJ =o,i 1 ... 1 n > 0, (n = 0, 1, 2, ...) (3) 

where Si+j is set to zero for odd i + j. P(9,r) is unique 
if, in addition, X^fTI'^nM] -1 ' 72 ™ = 00 (Carleman's crite- 
rion pq]). This is satisfied if P(9, r) falls off exponentially 
or faster as \9\ — > oo. 

Direct calculation of the moments of P(9, r) verifies 
that the r dependence implied by (0) is 



'■<».'>-f'(vK!"')* 



(4) 



where r and r' are both in the scaling range and p(a, R) 
has moments a 2n p(a, R)da = BS 2n . Explicitly, 



T(x,i) = K V 2 T(x,i). 



(1) 



p(a,R) 



4>(k,R) 



cos(fca)0(fc, R)dk, 



/(2n) 



(5) 



(6) 
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Eqs. (0)-@ are easily generalized to the case where the 
S2n(r)/S2n(r') do not scale as powers of rjr'. 

The Holder relations are equalities for normal scal- 
ing C2n = nC, 2 , which yields p(a 7 R) — 5 (a — R&/ 2 ). 
The smallest possible values of the higher C,2 n plausibly 
are C,2 n = C2, which yields p(a,R) = (1 — R^ 2 )5(a — 
(+0)) + R^ 2 S(a — 1). This corresponds to regions of 
constant T interspersed with regions where 9(r) is as 
large as the macroscale differences 9(L); ([TJ) cannot mag- 
nify differences in the injected T, only change their spa- 
tial scale and relax them. Unless the scaling is normal, 
p(a,R) increases in normalized width ajR^I 2 as R de- 
creases, so that P{6, r) cannot have a similarity form like 
r-PF{6/rP). 

In both cases Can — nCa and (2™ = C2, the tail oiP(9, r) 
is Gaussian if P(9, r) is Gaussian for r ~ L. In a Gaussian 
pdf lnS^n ~ nlnn for large n, to within terms of order 
n. If lnS^nM ~ nlnn for r ~ L, it follows that, for all 
"C2 > C2n > C2, lnS , 2„(r) ~ nlnn for large n throughout 
the scaling range. This is asymptotically Gaussian be- 
havior of the tail. An exponential or stretched exponen- 
tial pdf with tail In P(0,r) cx -\9\ b , with b < 2, instead 
requires lnS^n ~ (2n/6)lnn as n — * 00. Thus the tail 
of P(9, r) cannot fall off more slowly than Gaussian in 
the scaling range. This does not preclude exponential or 
slower fall-off in the dissipation range of r. The same re- 
sult holds for the pdf of velocity differences in the inertial 
range of Navier-Stokes turbulence. Other arguments for 
faster-than exponential tails in the inertial range have re- 
cently been given by Ching and Procaccia and by Noulez 
& 

In interpreting this result, a distinction must be made 
between the presence of a central cusp in the pdf, ex- 
pressing the existence of regions in which there is almost 
no excitation, and the shape of the far tail. Confusion 
between these two features can result in spurious fitting 
of a pdf to an exponential or stretched exponential. An 
example of a cusped pdf is that implied by the linear 
ansatz and shown in Fig. 2. 

If there is a range of r where the source term may 
be neglected, the steady-state balance equation for the 
even- integer structure functions S l 2„(r) is pj 



2 d 



Or 



„d-l 



.d-l 



n(r) 



dS 2n (r) 
dr 



nJ2n(r), 



(7) 



in the limit of rapidly-changing (white) velocity field. 
The left side of ^ is derived perturbatively from the 
u • V term in (|l|) in the white limit. Here d is space 
dimensionality, n(r) is the two-particle eddy-diffusivity 
scalar exerted by the velocity field, and 



J2 n {r) = 2n([9{v)] 2n - 1 H[9{v)}) 



H[6(r)]=((V 2 x + V 2 x ,)6(r)\6(r)), 



(8) 



(9) 



where (•|0(r)) denotes ensemble average conditioned on 
a given value 9(r), and x' = x + r. 

Power-law dependence of S2 (r) in the inertial range 
is assured from the exact, closed equation of motion for 
S2{r), provided that n(r) has the form 



fj(r) = ri(L)(r/L)* (0 < £ < 2), 



(10) 



for L ^> r 3> id- The steady-state scaling exponent for 
S 2 {r) found from (0) is (2 = 2 - £. 

There is no a priori assurance of power-law inertial- 
range scaling of 5 l 2„(r) for n > 1. If it exists, the balance 
of advection and n terms in (Q) implies that the J2 n {f) 
have the form 



Jin{r) = nC 2n J2(r)S2n(r)/ S 2 (r), 
where the C2n are dimensionless constants and 

(2n{(2n + d ~ ( 2 ) 



C 2 n — 



nd( 2 



The Holder inequalities (2™ < nC,2 give 

1 



C 2 n < 1 



<2 



(11) 



(12) 



(13) 



Equality in ( |l3| ) correponds to normal scaling (2™ = nC,2- 
The simplest candidate for an anomalous scaling solu- 
tion to (0) is the "linear ansatz" (!]]: 

H(9,r) = J 2 (r)6/S 2 (r). (14) 

It corresponds to G% n — 1 for all n and yields 

C2„ = yind(2 + (d- C2) 2 - \{d- C2). (15) 

Linear relations like (jl4|) were earlier proposed by Ching 
and Pope 

Using Mathematica® 3.0, I have verified numerically, 
for a number of values of d and (2 , that (|l5|) satisfies (|^) 
up to 2n = 100, even in the extreme case where P(9,r') 
in (|J) is taken as a 5 function. The C,2n corresponding to 
C211 = l + /3(n— 1) appear to satisfy (||) for < (3 < (2/d. 
This case [23 ] yields C,2 n c< n as n — > 00; (3 — (aid is nor- 
mal scaling. The perturbation analyses Q 0Ej-[L3 S^ ve 
nC,2 — Can that are quadratic in n for moderate n. They 
also satisfy (||); the full set of fan, and hence P(9, r), are 
not predicted. 

A conditional mean related to dissipation may be de- 
fined by 



K(9,r) = (\V X 9 + V X ,9\ 2 \9) 



(16) 



If r >■ id, the VzVz' crossterm in ( |l6| ) is negligible, and 
K{9,r) w (IV^^I 2 + \V X >6\ 2 \6), the conditional mean of 
the average of the dissipation at the points x and x'. 
H(9,r), K(9,r), and P(9,r) are related by an identity 
that follows from homogeneity alone pl|,p2[: 



2 



H(9,r)P(9,r) = — [K(9,r)P(9,r) 



(17) 



If the tail of P(8, r) is Gaussian, and K(9, r) is no 
stronger than a power of \9\, the leading term on the 
right of @ at large \8\ is oc -9K(9,r)P(9,r), whence 
H(0,r)/6<x-K(e,r). 

Monotonic growth of K(9,r) with \9\ at inertial-range 
r is a very weak, qualitative form of Kolmogorov refined 
similarity hypothesis If it holds, then the linear 

ansatz for H(9,r) and associated ( 2 n oc n 1 / 2 behavior 
at large n represent a lower bound for the asymptotic 
growth of C,2n- This bound is achieved if K(9, r) tends to 
a constant as |0| — » oo. Sublinear growth of H(9,r) at 
large \9\ would make K(9,r) a decreasing function of |0| 
and would yield slower growth of £ 2ri than n 1 / 2 at large 
n. Yakhot |l3| and Chertkov jl7j predict £ 2n — > constant 
asn-> oo, which can be shown to imply K(9, r) oc \9\~ 2 
as |0| -> oo [see @ and (||)]. 

The analytical relations between H(9,r), P(9,r), and 
the Q 2 n in the inertial range may be expressed in terms of 
a "co-pdf" P H {9,r) with moments ^ 0O 9 2n P H {9,r)d8 = 
C 2n S 2n (r) (C = 1). Then by (§) and"®, 



J7(0,r) _ J 2 (r) P H (0,r) 



9 



(18) 



5 a (r) P(0,r) 

and j°° oo P H {8,r)d8 = 1. If (|) is satisfied with S* 2 „(r) 
replaced by 5^(r) = C 2 „5 2 „(V}_, then P H and H(9,r)/9 
are everywhere positive. By (| 



^'^^c^^'^' 01 ' (19) 

Pn{9,r) is related to functions pjj(a,r) and <f>H{k,r) 
by equations identical to (Q)-(^|) except that a factor C 2 „ 
appears in ^>fj(k,r): 



b H (k,R) = l + J2 RC2 " C 2n(-k 2 ) n /(2ny.. 



(20) 



71=1 



Ph{9,t) = P(9,r) under the linear ansatz. Pjj(9,r) 
can be found analytically in some other cases that are 
formal solutions of (|l^). One is C 2 „ = 1 + (3{n — 1). The 
associated Ph is 



P H (9,r) 



1 



3/3 

T 



P(6»,r) 



/36> dP(0,r) 
"2 90 ' 



(21) 



Another case is £ 2 „ = C2 at all n, for which (^2|) yields 
C 2 „ = 1/n. These C 2 „ values are generated by 



P H (9,r) 



1 

W\ 



P(9',r)d9'. 



(22) 



|f| 



More generally, closed forms for Ph can be written if C 2ri 
is a finite polynomial in n. 

Eqs. (^|) and (|2^) converge rapidly because of the de- 
nominators (2n)L They have infinite radii of convergence 



if the ( 2 n satisfy the Holder inequalities. In some cases it 
is feasible to sum (||) and (^) numerically, perform the 
transforms and explore the shapes of P and Ph- Both 
4>{k,R) and <j>H{k,R) can be computed for large values 
of k, far out in the tails, by the use of Mathematica® 
to sum hundreds of terms at precisions of the order of a 
hundred decimal digits. 




FIG. 1. (j>(k,R) for d = 2, C2 = 1/2, R = 0.01 under the 
linear ansatz (solid line). The dashed line is <j)(k,R) if the 

linear-ansatz £an are reduced by n 1//10 , and the dotted line is 
4>H(k,R) for that case. 
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FIG. 2. P(9,r) for the linear-ansatz case (solid line); 9 rrns = 

10" 1/2 . The dashed curve is a normalized Gaussian. The 
inset is p(a, R) for the linear-ansatz case. 

Fig. 1 shows 4>(k,R) with R = 1/100 and £ 2 „ given 
by the linear-ansatz at d — 2, £ 2 = 1/2. Fig. 2 shows 
p(a, R) and P(0, r) as constructed from directly calcu- 
lated (f)(k, R) values for < k < 200 and an extrapolation 
for k > 200. The local exponent at k = 200 is -1.367 
and is slowly decreasing. P(9, r') is taken as a Gaussian. 

One can argue that, if the linear ansatz is self- 
consistent, the branch point in ( |l5| ) at 



c, 



-(d-C 2 )/2, q c = -(d-( 2 ) 2 /(2d( 2 ). (23) 
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should mark the limit of applicability of ([15]), and q c 
should represent the most-negative moment order that 
exists j35|. P(\9\,r') is here taken as a 5 function so that 
its moments of all negative orders exist. I therefore con- 
jecture that the ansatz gives (f>(k, R) cx fc~l 9c (k — > oo) so 
that p(a,R) cx a, I I 1 at the origin. The P(9,r) given by 
(|J) then is oc const. ~|6*|l ?c l _1 at the origin if 3 > \q c \ > 1 
and oc |6*|l 9c l _1 (infinite cusp) if |g c | < 1. The crossover is 
at = (2 — v3)<£ The data shown in Fig. 2 appear to 
be consistent with exponent —1.125 at k — oo, together 
with J Q <p(k, R)dk = 0, in accord with this conjecture. 

The primary fact here is that the linear ansatz ap- 
pears to yield a p(a, R) that is positive everywhere, en- 
suring positive P(9, r) also. The positivity of p(a, R) is a 
stronger result than (^). It implies positive P(9,r) even 
in the unphysical case where P(0,r') is a S function. In 
Fig. 2 there is an absolute cutoff of p(a) at finite a. This 
confirms that P(9, r) is Gaussian-like at large \0\, despite 
its appearance, if P(0, L) is Gaussian. 

Fig. 1, dashed curve is the 4>(k,R) that results from 
a gentle modification of these linear- ansatz exponents, 
C2n = Cfn" -1 ^ 10 ; where Cj^ n satisfies (|15|). This leaves £2 
unchanged and makes the C,2n go as nN^ as n — > 00. The 
large-fc (k ~ 200) tail of <j> appears to go approximately 
as fc-° 18 ; (§) is satisfied for 2n < 100. The implied be- 
havior of P at small 6 is 9~ Q1& . The result for 4>H:(k, r) is 
Fig. 1, dotted curve. Its large-fc behavior is qualitatively 
different from that of (f>(k, i?); 4>h continues to grow more 
negative at k = 200. Its tail is very different from that 
of cj>, which suggests bizarre behavior of H (9, r) at small 
9. The effective exponents C,§n = C271 + C<x n lni? for this 
case differ little from the ^n- If they are tried as expo- 
nents C211J they satisfy the Holder inequalities, but they 
violate @ at 2n = 56 if P{9, r') in (g) is taken as a S 
function. This illustrates that mild and innocent-looking 
deviations from the linear ansatz can destroy realizabil- 
ity. 

Numerical investigation is feasible for other values of 
£2 and R, and other deviations from the linear ansatz 
(e.g. C2n behaving as n 1 / 2 at large n, but with C2n 7^ !)• 
The results appear to be consistent with the conjecture 
that ( P3|) determines the behavior of P(9, r) at 9 = 
under (O). In particular, for d — 2, ( 2 = 1, R = 0.01, 
the largeTfc exponent of cf>(k, R) clearly lies between and 
— 1 (C^ — 0.536). Sensitivity to changes of the Cin away 
from unity seems to be widespread. 

I am grateful to Shiyi Chen, E. Ching, G. Eyink, T. 
Gotoh, I. Procaccia, K. Sreenivasan, and V. Yakhot for 
helpful interactions. This work was supported by The 
Department of Energy. 

* Internet address: rhk@lanl.gov 
[1] R. H. Kraichnan, Phys. Rev. Lett. 72, 1016 (1994). 
[2] V. S. L'vov, I. Procaccia, and A. L. Fairhall, Phys. Rev 

E, 50, 4684 (1994). 
[3] B. Shraiman and E. D. Siggia, Phys. Rev. E, 49, 2912 

(1994). 



[4; 
[5 

\-\ 

[§ 
[9 
[10 

[11 

[12 
[13 
[14 

[15 

[16 

[17 

— 

[19 
[20 



[21 

[22' 
[23 

[24 
[25' 
[26 
[27 
[28 



[29 

[30 
[31 

[32 
[33 



[34 
[35 



E. Balkovsky, M. Chertkov, I. Kolokolov, and V. Lebe- 
dev, JETP Lett. 61, 1012 (1995). 

M. Chertkov, G. Falkovich, I. Kolokolov, and V. Lebe- 
dev, Phys. Rev. E, 51, 5609 (1995). 

M. Chertkov, G. Falkovich, I. Kolokolov, and V. Lebe- 
dev, Phys. Rev. E, 52, 4924 (1995). 

M. Chertkov and G. Falkovich, Phys. Rev. Lett. 76, 2706 
(1996). 

A. Fairhall, O. Gat, V. L'vov, and I. Procaccia, Phys. 
Rev E, 53, 3518 (1996). 

K. Gawgdzki and A. Kupiainen, Phys. Rev. Lett. 75, 
3834 (1995). 

R. H. Kraichnan, V. Yakhot, and S. Chen, Phys. Rev. 
Lett. 75, 240 (1995). 

V. L'vov and I. Procaccia, Phys. Rev. Lett. 76, 2898 
(1996). 

A. Pumir, preprint 1995. 

V. Yakhot, submitted to Phys. Rev. E, 1995 and 1996. 
D. Bernard, K. G awedzki, and A. Kupiainen, preprint 



chao-dyn/960101S . 

M. Chertkov, G. Falkovich, and V. Lebedev, Phys. Rev 
Lett. 76, 3707 (1996). 

D. Segel, V. L'vov, and I. Procaccia, Phys. Rev. Lett. 76 
1828 (1996). 

M. Chertkov, submitted to Phys. Rev. E (:hao- 



Falkovich, I Kolokolov, and V. Lebedev, 



and 



dyn/960601l| ). 
E. Balkovsky, G. 

submitted to Phys. Rev. E ( |shao-dyn/9603015| ). 
E. S. C. Ching, Phys. Rev. E 53, 5899 (1996) 
E. S. C. Ching, V. S. L'vov, Evgeni Podivilov 
I. Procaccia Phys. Rev. E 54, 6364 (1996) ( |chao- 
dyn/9608008j ). 

O. Gat, V. S. L'vov, E. Podivilov, and I. Procaccia, sub- 
mitted to Phys. Rev. Lett. 
G. L. Eyink, Phys. Rev. E 54, 1497 (1996). 
E. S. C. Ching, V. S. L'vov, and I. Procaccia, Phys. Rev. 



E 54, R4520 (1996) ( |chao-dyn/9606017| ). 

G. Eyink and J. Xin, Phys. Rev. Le tt. 77, 2674 (1996) . 
R. Benzi, L. Biferale, and A. Wi rth (bhao-dvn/9702o"o5[ ). 
M. Vergassola and A. Mazzino ( |chao-dyn/9702014| ). 

S. Chen and R. H. Kraichnan, in preparation 

H. S. Wall, Analytic Theory of Continued Fractions (New 
York, Van Nostrand, 1948), Sees. 86-88. Example: the 
modified Gaussian moments 0.9, 3, 15, 105, 945, ... sat- 
isfy all Holder inequalities but violate (^) at 2n = 22. 
A. Noulez, unpublished; E. S. C. Ching and I. Procaccia, 
unpublished. 

E. S. C. Ching, Phys. Rev. Lett. 70, 283 (1993). 

S. B. Pope and E. S. C. Ching, Phys. Fluids A, 5, 1529 

(1993). 

E. S. C. Ching and R. H. Kraichnan, in preparation 
G. Stolovitzky and K. R. Sreenivasan, Rev. Mod. Phys. 
66, 229 (1994); G. Stolovitzky, P. Kailasnath, and K. R. 
Sreeivasan, J. Fluid Mech. 297, 275 (1995). 
Y. Zhu, R. A. Antonia, and I. Hosokawa, Phys. Fluids 7, 
1637 (1995). 

V. Yakhot, private communication. 



4 



